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Abstract 

Unidirectional start-up flow of a viscoelastic fluid in a pipe with frac- 
tional MaxweU's model is studied. The flow starting from rest is driven by 
a constant pressure gradient in an infinite long straight pipe. By employ- 
ing the method of variable separations and Heaviside operational calculus, 
we obtain the exact solution, from which the flow characteristics are in- 
vestigated. It is found that the start-up motion of fractional Maxwell's 
fluid with parameters a and /?, tends to be at rest as time goes to infinity, 
except the case of /3 = 1. This observation, which also can be predicted 
from the mechanics analogue of fractional Maxwell's model, agrees with 
the classical work of Friedrich and it indicates fractional Maxwell's fluid 
presents solid-like behavior if /3 7^ 1 and fluid-like behavior if jS — 1. For 
an arbitrary viscoelastic model, a conjecture is proposed to give an intu- 
itive way judging whether it presents fluid-like or solid-like behavior. Also 
oscillations may occur before the fluid tends to the asymptotic behavior 
stated above, which is a common phenomenon for viscoelastic fluids. 

Keywords: Viscoelastic fluid; fractional Maxwell's model; start-up flow; pipe flow; 
Heaviside operational calculus. 



1 Introduction 



'All things are movable and in a fluid state', which is a famous quotation from 
Thales of Miletos, the first philosopher of ancient Greece. 

Indeed, besides the most familiar fluids such as water and gas, most mate- 
rials in nature and industry, such as milk, oil, lava, etc., can be treated and 
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investigated as fluids. However, some of them not only have the viscosity like 
Newtonian fluid, but also exhibit Hooke's elasticity. They are viscoelastic ma- 
terials. Different models of viscoelastic materials were obtained and studied 
during the past hundreds of years; for example. Maxwell's model, constructed 
by a spring and a dashpot in serial, was largely investigated in the last century. 

Recently, fractional Maxwell's model, constructed by two fractional element 
models in serial, attracts a lot of researchers' interests [I] [2] [3] [4] [5] [15] [18] [19] [20] . 
Let (7 be the shear stress and e be the shear strain. The constitutive equation 
for fractional Maxwell's model is given by 

fT + A"— c7 = ^A^— .e, 0<a</3<l, (1) 
dtp ' - -H - , \ J 

where E the shear modulus, and A is the relaxation time0 
In the case of a = 0, equation ([Ij degenerates to 

which is just the constitutive relation of a fractional element model with the 
shear modulus E jl. In fact, the usual expression for a fractional element model 
was first introduced by Scott Blair [6] [7] 

. = i^.A^^e, (2) 

where E^ is the shear modulus. The mechanics analogue of a fractional element 
model can be found in [T][5][H][Tn]. 

In the case of a = /3 = 1, equation (|T]) degenerates to the constitutive 
relationship of classical Maxwell's model. 

Physically, fractional Maxwell's model can be considered as two fractional 
element models in serial, with orders 71 and 72 satisfying 

a = l7i - 72I, = max{ 71,72}. (3) 

Fig.l gives the mechanics analogue of fractional Maxwell's model, where a tri- 
angle denotes a fractional element model. 

The continuous interest of fractional Maxwell's model is perhaps due to the 
special known and unknown property of this model and due to the rapid de- 
velopment of fractional calculus during the last fifty years. Palade et. al. in 
[15] derived fractional Maxwell's model from the linearization of the objective 
equation and discovered the anomalous stability behavior of the rest state in 
three dimensions. Tan et. al. in [3 [9 studied four unsteady flows of a vis- 
coelastic fluid with generalized Maxwell's model between two infinite parallel 
plates. Vieru et. al. in [T5] studied flow of a generalized Oldroyd-B fluid due 
to a constantly accelerating plate, which includes generalized Maxwell's model 

^Some studies on fractional Maxwell's model only concern the particular case with /3 = 1, 
i.e., the so called generalized Maxwell's model. 
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Figure 1: Mechanics analogue of fractional Maxwell's model. A triangle denotes 
a fractional element model. 

as a limiting case. Hayat et al. f5] also discussed fractional Maxwell's model 
and studied three types of unidirectional flows which were induced by general 
periodic oscillations of a plate. Hernandez- Jimenez et. al. gave some experi- 
mental results for oscillating flows with fractional Maxwell's model '4^, which 
encourage more studies. Yin and Zhu [1] studied the oscillating flow with frac- 
tional Maxwell's model in an infinitely long pipe and found interesting results, 
for example, the resonance peaks was discovered to be different with those of 
ordinary Maxwell's model. 

A lot of interests and studies were also given to the unidirectional start-up 
pipe flows, which has a significant practical and mathematical meaning. Zhu and 
Lu et. al. in [M] studied characteristics of the velocity filed and the shear stress 
field for an ordinary Maxwell 's fluid and discovered the oscillation phenomenon. 
Fetecau in [171 studied the analytic solution for an ordinary Oldroyd-B fluid and 
several limiting cases such as ordinary Maxwell's fluid; the velocity profiles for 
the steady state are the same in all types of fluid they studied. Further, Tong 
et. al. in jl6 l studied the exact solution for the fractional Oldroyd-B model in 
an annular pipe by using Hankel-Laplace transform. Using similar but different 
methods, Zhu and Yang et. al. in [13] studied the exact solution and flow 
characteristics for the fractional element mode^ with parameter /?, the most 
fundamental model in all fractional derivative models. They found oscillation 
phenomenon and solid-like behavior for certain /3. 

As far as we know, the characteristics of start-up pipe flow with fractional 
Maxwell's model have not been well-studied yet. In this paper, we investigate 
basic characteristics of such flows through studying the exact solution. The 
fluid is quiescent in the beginning in an infinitely long pipe, and then it will 
be suddenly started by a pressure gradient which remains constant after the 
starting moment. 

The start-up pipe flow with a dashpot model (i.e. Newtonian fluid) is clas- 
sical in fluid dynamics; with ordinary Maxwell's model the motion would be 

■^We mention that the fractional Oldroyd-B model in the study of Tong et. al. does not 
include the fractional element model as a limiting case. 
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very interesting because of the occurrence of oscillations [13]; and with the 
fractional element model solid-like behavior was discovered [13] as we already 
mentioned above. One must be curious what the interesting phenomena for 
fractional Maxwell's model are and whether we can summarize a way to deduce 
the characteristic of a certain viscoelastic model if just given the constitutive 
relation. 

This paper is organized as follows. In section 2, we present the governing 
equations and initial boundary conditions, and solve this initial boundary value 
problem through the method of variable separation and Heaviside's operational 
calculus. In section 3, we discuss the flow characteristics and give some obser- 
vations. In section 4, we make the conclusion and propose a conjecture. 



2 Governing equation and exact solution 

We consider the start-up flow of a fractional Maxwell fluid in an infinitely long 
pipe with the radius a. In the beginning, the fluid in the pipe is at rest; then 
it is suddenly started by a constant pressure gradient. Take the pipe axial 
direction as the z-axis. We construct a column coordinate system (r, 0, z) and 
let V{r,9,z,t) denote the flow field. 

Since the flow is axisymmetric, we assume V only has the axial component 
and does not depend on 9, i.e., 

V = uir,t)e,. (4) 
The governing equations of the motion are given by the continuous equation 

V • y = 0, (5) 
as well as the momentum equation 

dV 

p— = -Vp + V • a, (6) 

where ^ is the material derivative, p is the density of the fluid, p is the pressure 
fleld and a is the stress tensor field. 

The start-up flow considered is driven by the pressure gradient field given 

by 

^(t) = -Gh{t), (7) 
oz 

where G is a constant and h{t) is the Heaviside function, defined by h(t) = 
0, t<0; h{t) = 1, t > 0. 

The constitutive equation (|T]) gives: 

The initial and boundary conditions are given by 

u{a,t)^0, u(0,i)<oo, u(r,0) — (r,0) = 0, t > 0, < r < a. (9) 
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Substituting ^ to equations ^ and ([S]), and considering ([7]) and (H)), we 
find 

pfi + # = fi + a"|^)g/. + i^A^m f + fj). (10) 

'^V dt"') dt V at"/ dtP-^\rdr dr^ J ^ ' 

Let the radius a be the characteristic length, pa^ /E\ be the characteristic 
time, and Ga^ /EX be the characteristic velocity. By simple algebraic manip- 
ulations, we get the dimensionless governing equation and the dimensionless 
initial-boundary conditiontlfl: 



< r < 1, i > 0, 

and 

w(l,t) = 0, u(0,i)<cx3, t>0, (12a) 

w(r,0) = — (r,t)=0, 0<r<l. (12b) 

We use the method of variable separation and Heaviside operational calculus 
solving equations ([TT|) and Let 

M(r,i) = w(r)T(t). (13) 

Substituting this expression to the homogenous equation of equation (jlip . we 
obtain 

-(-.'(,•)+."(.)) . \;._,^,,,(,/ ^ - (14) 

where k is some appropriate constant to be determined. Solving the eigenvalue 
problem: 

-v'{r)+v"(r) + k^v^Q, (15a) 
r 

v{\) ^ 0, v(0) < oo, (15b) 
we obtain the discrete eigenvalues 

fci < fc2 < fca < (16) 
as well as the corresponding eigenfunctions: 

Vm^ Joikrar), TO =1,2,..., (17) 



•^We still use notations u,r,t but here they denote dimensionless quantities. In the follow- 
ings we only consider dimensionless quantities, so this change of notations will not bring any 
confusions. 
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where km is the rrith positive root of the zeroth Bessel function. 
The final solution is constructed by 



u{r,t) = J2 AmTm{t)Mk^r), (18) 

m— 1 

where m — 1,2, 3,..., are constants to be determined and T^it) are func- 
tions of t to be determined. Substituting this expression of u{r,t) to equation 
(fTT|) . we find 

£ Am{KM+>^''T"^\t)+kiy-'T!:-\t))Mkmr) = hit) + ^p-^; (19) 

m— 1 ^ ^ 



Since 



Ht) + 4^ - (hit) + V 

Til -a) V ^ ^ ril-a)J ^^hnJi 



ikm) 



(0 < r < 1), (20) 



comparing the coefficients of the eigenfunctions appearing in equations (|19[) and 
(1201) we have 

as well as 

r^(t) + A"T„"+i(t) + fc2,A^-iT,fr^(i) = hit) + ^—-y (22) 

r(l - a) 

We solve equation (|22p by applying Heaviside operational calculufQ. Let 
p — and let Tmit) ~ Yhit) where Y is an operator to be determined. Noting 
that 

vMt) . (23) 



we have 



As a result, 



Y = ^ + (24) 



1 + A"p" 
p + A"p'^+i + k^^Xi^'^p' 

By the definition of the Heaviside operator [ll , it yields 



Trnit) = ^i, Mt). (25) 



''Heaviside operational calculus is in fact equivalent to Laplace transform method, but the 
former method is more intuitive: The spectral parameter has a clear meaning. 
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where L is a contour in complex z-plane parallel with the imaginary axis, and 
is determined by the requirement that there be no singularities of the integrant 
on the right of L. 

The final solution we construct is given by 



uir t) ^— V '^"^^-'^^ / ^ + ^"^^ e-dz (27) 



Substituting this equation to equation (|5]) and assuming the shear stress field is 
zero at i = 0, we obtain 

_ EXP ^ -2Ji(fc^r) f , 

^"""2^^^^, Ji(fc„) Az2 + A"z"+2 + fc2 A^-i^/^e dz. 

In the case of a = 0, i.e., the case of the fractional element model, solution 
(P7|) reduces to 

u(r,t] ^ = > - — ^-—T / — , „ ,o 1 a c'^^dz. (29) 

^ ^ 2^V^ kmJiikm) Jl + klX^-^zP ^ ^ 

By an inverse formula given in |10|[p.d. 271-273], we can simplify expression 
HMD: 

"(^'*) = E |^^^^ti?2-,,2(-fc,lA^-4^-^/2), (30) 

m — l ^ ^ 

where £2-/3,2 is the Mittag-Leffler function. Note that the Newtonian fluid 
corresponds to the particular case of (a = 0, /3 = 1). Substituting (3 — 1 to 
equation (1301) . we obtain 

(31) 

which is the classical dimensionless solution. 

In the case of a = /3 = 1, i.e., the case of ordinary Maxwell's model, our 
solution agrees with the solution obtained by Zhu et. al. |14) . 



3 Results and discussions 

Due to the simplicity of the fractional element model, which also plays the 
fundamental role of constructing different fractional models, we first recall some 
results in [13] of start-up pipe fiow for the fractional element model (Scott Blair's 
model). Without lost of generality, we take A = 1 in equation (pO)) . As Fig. 2 
shows, for the case 

< ^ < 1, 
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oscillations occur just after the fluid is started; the smaller the parameter j3 is, 
the stronger the elasticity is. And as t goes to infinity, the center velocity tends 
to be 0. In fact, recall the asymptotic formula for £^2-^,2 [lOj : 



1 1 



r(/3)z 



+00, 



(32) 



We know from this formula that the series expression pO[) is uniformly conver- 
gent at least for any fixed r. Let t — > 00; we obtain that 



lim u{r, t) = 0, 



< r < 1. 



(33) 



The only exception is the case of /3 = 1, i.e. the case of Newtonian fluid. No 
oscillations would occur and as t goes to infinity, the center velocity will tend 
to be a steady constant 0.5. In this classical case. 



1 



lim u(r,t) = -(1 - r-^), < r < 1. 



(34) 



Fig. 3a) and 3b) shows the velocity profiles of different t, with /3 — 0.4 and 




Figure 2: Center velocity with respect to t for the fractional element model. 

(3 — 1 respectively, which gives intuitive pictures in mind. 

Based on these discussions of of the fractional element model, the mechanics 
analogue of fractional Maxwell's model (see Fig.l) would help us do further 
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a. /? = 0.4 



h. 13 = 1 



Figure 3: Velocity profiles at different t for the fractional element model. 

investigates. This idea will be directly applied in the following discussions. Still 
we will take A = 1 without loss of generality. And we discuss the start-up flows 
with different a and /3 through studying equation (j27l) . 

Fig. 4 gives the curves of the center velocity with respect to t when a — 0.6. It 
can be seen that for 0.6 < /? < 1, the center velocity of the pipe will tend to be 0, 
which means the corresponding fractional Maxwell's model will finally represents 
solid-like property; for (3 = 1, the center velocity will tend a constant 0.25 and 
the fractional Maxwell model finally represents fluid-like property. Furthermore, 
the smaller /3 is, i.e., the stronger the elasticity of the larger order fractional 
element is, the stronger the oscillating phenomenon is. 

Fig. 5 is similar with Fig. 4 but the phenomenon is more clear. It gives the 
relation curve of center velocity with respect to t with different a, when /? = 0.6. 
It can be seen that, as t goes to infinity, the center velocity of the pipe tends to 
be 0, which indicates the corresponding fractional Maxwell's model represents 
solid-like property. When a increases, i.e. the difference between the orders of 
the two fractional elements increases, or say the elasticity of the smaller order 
fractional element strengthens, the oscillating phenomenon becomes stronger, 
which helps the fluid to be at rest, although in short time it accelerates the fluid 
more. 

Fig. 6 shows the relation curve of the center velocity with respect to t in the 
case of /? = 1. We see that if < a < 1, the center velocity will tend to a 
constant 0.25 as t goes to infinity; if a = 0, i.e., the case of Newtonian fluid, the 
center velocity will tend to 0.5. In both cases, fractional Maxwell's fluid will 
represent fluid-like property for large t. And when a 0, oscillation also occurs 
which means the fluid also has elasticity, however, since /? = 1, the fluid-like 
property will lead the way. 

It deserves to point out that these results partly agree with Friedrich's [12] , 
who first pointed out that the fractional Maxwell model represents solid-like 
property as long as /3 < 1. 

An intuitive way to deduce whether a fractional model would exhibit solid- 
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Figure 4: Center velocity with respect to t for the fractional Maxwell model in 
the case of a = 0.6. 



like or fluid-like property is to consider the fractance of the mechanics analogue 
of the model. According to the analysis of the fractal construction of these 
models [IHlIll, we can indeed deduce solid- like property without calculations: A 
fractional model (which itself is not a string) presents solid-like behavior if and 
only if there exists a spring path from one side to the other in any spring-dashpot 
fractance of its mechanics analogues. 



4 Conclusions 

An exact solution of the start-up pipe flow with fractional Maxwell's model is 
obtained and the flow characteristics are discussed. 

In the case of /3 ^ 1, the motion of fractional Maxwell's fluid in a pipe tends 
to be at rest as time goes to infinity; otherwise if /? = 1 the flow will have a 
parabolic-like profile as t goes to oo. 

Indeed, for a Scott-Blair's model (fractional element model) with parameter 
(3, as long as /3 < 1, the model presents a solid- like property. This was showed 
by Zhu and Yang et. al. in p3] in the study of start-up pipe flow. 

For a fractional Maxwell model, i.e., two fractional element models in a 
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serial connection, the fluid will present a solid-like property as long as the two 
fractional elements both present solid-like behaviors. On the other hand, if 
one of the fractional elements in serial is Newtonian, even the other one is a 
spring, the serial compound will be of fluid-like property and the flow will keep 
a stationary velocity profile after infinitely long time. Our result agrees with 
the results in [T^], but from different points of views. 

Moreover, the stronger the elasticity of any of the fractional element models 
in the serial connection is, the stronger the oscillation of the corresponding 
fractional Maxwell model is. 

From these results for the case of fractional Maxwell's model, we conjecture 
that a viscoelastic model presents a solid-like behavior if for any spring-dashpot 
fractance of its mechanics analogues, there exists a spring •path(a, path with only 
springs) from one side to the other. This conjecture, if proved, will generalize 
Friedrich's result in [T^ . 

Acknowledgements We are very grateful to the referee for suggestions and 
constructive comments. 



11 




Figure 6: Center velocity with respect to t for the fractional Maxwell model in 
the case of /3 = 1. 
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